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The dominance of xo in the data for e'*'e~ ip + X/(cc) is shown to violate OZI factoriza- 
tion. Single gluon exchange gives a non-factorizing effective S ■ L interaction that generates a large 
scalar production amplitude. This also has observable effects near threshold in e'^e~ — > D^*^Dj, 
where enhancements of D* Do and DDi channels are predicted. Further tests and implications are 
discussed. 



I. INTRODUCTION 



Hadron decays occur dominantly when the colour sources (such as QQ) are separated by 0{Aqcd) such that the 
energy in the fields is enough to enable the creation of a light qq pair. This has been embodied in the phenomenological 
OZI rule [H, incorporated in models and more recently studied in lattice QCD Q . It has been shown Q that lattice 
QCD appears to confirm assumptions implicit in many models of strong decays, namely that there is a factorization 
of the constituent spin S and total J of the hadrons (in models this equates to a factorization of S and L) . Similar 
conclusions and results are being found in the AdS/QCD correspondence @. One consequence of factorization is that 
e+e" [0++] must be smaller than at least one of e+e" and e+e" [2++] [J. This is in 

marked contrast to data on e+e~ ^ tp + X where the xo appears to dominate as reported in ref 

Decays that involve the creation of heavy flavours where 2mQ > Aqcd (e.g. cc —^ cc + cc) may differ radically from 
those involving creation of light flavours. The OZI process is suppressed because 2mQ > Aqcd would require the 
colour fields of force to extend over excessive distances without having created light qq. This is highly improbable. 
A way for such decays to be triggered is if the required energy is supplied by a hard process such as single gluon 
exchange. As the gluon will in general transmit information about S and L from the initial quarks to the created 
pair, we expect that factorization will not occur in such processes. 

Braaten and Lee 7] present a calculation of e+e^ ''i^ + XJ which breaks the limits imposed by factorization, 
and which appears to be supported by the data - in particular the xo dominates the xi and X2- This situation has 
drawn little comment either theoretically or experimentally despite its important ramifications for strong interaction 
phenomenology. The purpose of our paper is to: draw attention to the importance of Braaten and Lee's results, 
expose the origin of their remarkable result given their assumptions, extend their result to other cases and develop 
stringent experimental tests of the OgE mechanism as opposed to factorized interactions. 

In section 2 we review the predictions of factorization. In section 3 we study the general structure of the gluon- 
driven processes producing -0 + xj as a function of J and expose the origin of the large xo amplitude in ref. 7]. We 
apply these results to the production of pairs of charmed mesons and find that they lead to an enhancement of DDi 
and D*Dq channels near threshold, but make no contribution to D*Di or D*D2- Further, we predict that the helicity 
amplitude I)*(±1)D2(=f2) = 0, whereas ';/'(il)X2(T2) is the largest of all the ■0X2 amplitudes. 
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II. FACTORIZATION 



A feature common to many models of hadron decays is that the constituent spins factorize from the total angular 
momentum of the hadrons. Ref [1] showed that this feature appears to be confirmed by lattice QCD and refs [3| 
showed how this property explains relations among various amplitudes that are common to many many specific 
models. This factorization property (eq. 5 of ref underpins the empirical result that hadron loops (such as in 
cc DD,DD*,D*D,D*D* — > cc) give universal mass shifts throughout an L multiplet such as: ft-c,X0jXijX2 @]- 
It is also claimed to be consistent with AdS/QCD correspondence^ 

Ref Q derived various consequences of factorization for strong decays of hadrons. In particular this work showed 
that factorization constrains the relative populations of final states in e+e~ — > V^i^S*!] + {S, A,T)['^Pj]. 

In the factorization scheme, the decay of a transversely polarized ^Si —^^ Si +^ P2, with the tensor meson maximally 
polarized along the decay axis is predicted to vanish (see eq. 34 in Q). This selection rule is a particular test of 
factorization. Also ref Q found that for ^5*1 decay, the relative rates in S'-wave, 5^, and I?-wave, are: 



^Si+^Po : ^Si+^Pi : ^Si+3P2 : ^Sq +^ Pi = : iS^ + : 6D^ : S^ + D^. (1) 

Among various consequences, of interest to the present discussion is that 

a(3Si ^3 Si +3 Pi) = la(3Si ^3 s, +3 p^) + l^^s^ +3 p^) (2) 

3 6 

and hence 

a{e+e- Si Si +^ Pi) > (j{e+e- -^^ Si Si +^ Pq). (3) 

One of the central applications of the present paper will be to test these predictions against data on e~^e~ — ^ '0 + XJ 
where preliminary indications are that the relation eq. [3] is violated f6| . 

The constraints of factorization become more powerful near threshold where S'-wave dominates. For a 3Si initial 
state 

/3r 



<^rsi~*^x2) (4) 
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^fSi-^^Xo) = ^fT(3Si^Vxi) 



Analogously, for a 3Di initial state 



a(3Di ^ 0xo) ^ (5) 







a(3Di^0Xi) - |(7(3Di^te). 

Finally one may allow for a coherent mixture of 3Si and 3Di initial state. Results become model dependent but 
cr{'ipxo) cannot be made larger than both criipxi) cr(V'X2)- Thus in the region of threshold factorization forbids a 
dominant (t(V'Xo)- 

This is interesting in view of the data on e+e~ — s- + A" at 10.6 GeV cm. energy, which show three prominent 
enhancements X in e'^e" — > ip + X Q, consistent with being the rjc, and xo- The observed pattern of states appears 
radically different to what is seen for light flavours: the apparent prominence of Yn with only a hint of xi and much 
suppressed xi contrasts with light flavours where e+e" — s- w/2 is clearly seen [10l |. The charmonium data ^ are 
significantly above threshold and so our general restrictions against e^e~ ^ + Xo need not apply. However, we 
shall see that the large rate for this channel is a signal for factorization breakdown even away from threshold, and 
inspires the question: what is required to create a dominant e~ — >■ V' + Xo amplitude? 



III. GLUON EXCHANGE STRUCTURE 



A. Non-relativistic reduction of the Feynman Amplitude 



If we require in fig. 1 that the upper cc pair produce the C = — meson then at leading order four diagrams 
contribute to e~'"e~ — > t\}XJ- These consist of gluon emission from c or c in either of two topologies: exchange within 
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the tp (figs. 1(a), 1(c)) or within the x (figs. 1(b), 1(d)). To contrast with the factorization amphtudes most directly, 
we shall restrict our attention to the threshold region. As in ref. we set the masses of the tp and x each = 2mc- 

We first make a non-relativistic reduction of the Feynman amplitudes into 2x2 block matrices sandwiched between 
two-component spinors. In so doing it is important to note the role of the (anti-)fermion propagator between the 
virtual photon and the gluon. 

Matrix elements in the explicit non-relativistic limit are discussed in a consistent phase convention by Ackleh et 
al. which we adopt here (see Appendix B of ref |ll|, especially eqs. B5-B7). Care is required to track phases 
and so we define here our choice of some convention dependent quantities. Dirac spinors are normalized to 1; 
J e~*^P~'''''^d^a; = {2tt)'^5'^{p — q); h = c = 1; four vectors are denoted by, say, p while three vectors are given by p; 
our metric is (-1-,-,-,-). Particle Pauli spinors, ^, and anti-particle Pauli spinors, rj are explicitly: 

i+=V- (6) 



The momentum routing is defined in fig. 1. 

The analysis of Ackleh et al. was for a qq wavefunction of well defined "^^^^Lj with no constituent propagator 
effects considered. In e^e~ annihilation of the present paper, we explicitly consider the constituent propagator 
between the photon and exchanged gluon. This has the effect of introducing more than just a sing le ^■^+^Lj state for 
the photo-produced cc, i.e. it is not simply ■^Si. 

From the usual Feynman rules, the matrix element for fig. 1(a) is: 

= uiigs-/"—^ "2 —teec-f''v4-Y^U3igsj^'v2. 7) 

where gs is the strong coupling constant and Ec is the ratio of the electron and c quark charges and e is the electron 
charge. Our primary interest is in the J-dependence of the production of ^Si -|-^ Pj near threshold. 

In NRQCD the appropriate L state is extracted by expanding the matrix clement into a double Taylor series of 
the internal momenta qi and q2 (about 0) and selecting the appropriate power of the internal momenta. The angular 
dependence of the appropriate term in the Taylor series is recombined with the spin of the hadron using covariant 
generalizations of Clebsh-Gordon coeffecients and the radial dependence is absorbed into a model-dependent vacuum 
saturated analogues of the NRQCD matrix element, (Oi),/, say. The (Oi)3pj are approximately independent of J, as 
one would expect from any spherically symmetric hadronization process, and as such are collectively denoted (Oi)3p. 

For the S and P-wave mesons of interest, the appropriate term is the one with a single power of q2 and zero powers 
of qi which is ensured by setting qi — 0. Since we are interested in the threshold limit, we simplify the series by 
taking Pi — P2 = 0. After expansion and algebraic simplification of the matrix element, the q2 will be replaced with 
and its radial dependence absorbed into the vacuum saturated matrix element. Finally the spin of the ^Pj and 
its orbital angular momentum will be combined to obtain the appropriate J. This is the non-relativistic reduction of 
the covariant NRQCD technique for ■^Si -l-'^ Pj. 

The matrix element is effectively simplified to: 
M 




where 92 = Q2/'2mc 

We generalise the above to allow the quarks of the flavour produced by the gluon to have mass m f . Then we define 
6 = 2mf /rric] in the particular case of ■0 -I- xj, = "^c and hence 5 = 2. For flavoured mesons with my — s- one 
has (5^0. We shall retain the general form so that applications to other combinations of flavours may be made. The 
q2 = <72/2toc then generalises trivially to contain the relevant my in the denominator: (72/ = g2/2my = |g2. N also 
has an implicit dependence on 5 but the dependence does not matter for the relative rates we are interested in. Hence 



LY7V('?')^2 (8) 



92/ 
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We can split the equation into a time-like component (7''7°) and a space-like component (77). Taking the leading 
order terms for the time- like component, 



Having written the matrix element in terms of block 2x2 matrices and particle, anti-particle Pauli spinors 
we can use 

(1, A; i, A) ^ -^^-e* {S) (10) 

whereby the substitution, r/^^ __i=o-.£* projects the matrix element into the spin triplet state. The terms of 

exactly one power of 92 are retained {qi — ensures that only terms with zero powers of qi were kept). 92 is replaced 
by el and its radial dependence absorbed into the vacuum saturated analogues of the NRQCD matrix elements. The 
result is the projection of the matrix element into the ^Si +^ P state: 

and then using 

iTr {aV^ff'^cr'} = 6'^ 6''^ - 6'''S^' + 6'^5^^ (11) 

we obtain 

^3 ^3 p) ^ n'^-^{0,U,{0^).^ ( el.ael.el - el-ae^-el + el-ae^l ) . (12) 

where a is the polarisation of the photon, n is the polarisation of the spatial components of the gluon, is the orbital 
angular momentum of the "^P meson and ei,e2 ^^^e spin polarisation tensors of the ■^Si,'^P mesons respectively. 
The time-like components of the amplitudes cancel among the various graphs. Therefore we do not consider them 
further and focus on the space-like components. 

For the space-like component, we have 



= _7VTr 



l\ ( cr"\ (mc{2 + 5) -(T.q\ f a'' 
q) i^+cr" ) \ cr.q -mj) \-a'' 

M^^-'ir Si P) = -i7V(Oi)3s, (Oi)3pTr | (a. ela.a - ^-n.aa.el ) rT.e\rT.e\ | . (13) 
Then as before, using ea llli this may be written 

M-^^(^* ^3 ^3 p) ^ Ar(o^)3s^(Oi)3p i ^-±le\.ael.e^ ^el.ae^.el - ^e^.ae^.ei ] . (14) 



The amplitude in fig. 1(b) is calculated in similar fashion and gives 



(15) 



2 

Explicit calculation of the other diagrams confirms the symmetries 
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Ma7«7''(^* ^3 ^3 p) ^ ^C7°7''(^* ^3 ^3 p) 

= -X'^^V(7* ^3 ^3 p) ^ _^d7°7''(^* ^3 ^3 p) 

>l^-y^(7* ->3 Si +3 P) ^ M'^'^'^i-f* Si +3 P); 

^3 ^3 p) ^ yV^dT-T-^^, ^3 g^ _^3 p) (^g) 

Thus the total ampUtude at threshold becomes 

Mir Si +3 P) ^ N{Oi)3s, (Oi)3p f (25 + 4)elael.e; - ^el-ae^.e^ ) (17) 



For arbitrary final state meson momentum the amplitude is 

M{j*" Si +3 P) cx Ami £1 ■ £2{6PiPi ■ £l ~ sel) 

+ £2 • eLs(s£i - 2Ff £1 • P2) + 24ml Pi ' £l{£2£i ' ^2 - £^£2 • Pi)- (18) 

a is the photon's helicity. We are grateful to Dr Jungil Lee for providing this expression. As in the threshold limit, 
this indeed reduces to our eq. [17] with 5 = 2 for physical values of the photon helicity a. 

Of the two contributions in eq. 1171 only the e*j^.ae\.e2 factorizes L and S. The first term couples Lip ^ and Ss^ j 
and survives everywhere except the unphysical case of threshold for 5 = —2. In the physical cases of interest, this 
term turns out to dominate. 



B. Amplitudes in NRQCD 



In the particular limit of &t threshold, 5 = 2 the amplitude is: 

M(r + cc(3p)) = 27V(Oi)3s^ (Oi)3p ( 4el.a£l.£*^ - el.ael.el ) . (19) 

For charm pair production at threshold where (5 the amplitude simplifies to: 

Mij* ^D*+ cqCP)) = N{Oi)ss, (Oi)3p ( 4£*.aei.e; ) . (20) 

Note that N and the vacuum saturated matrix elements will be different in the charm pair case as opposed to the 
double charmonium case, but this will not matter for relative amplitudes in the D* Dj channels of interest. 

Eq. [18] is the starting point from which the amplitudes of ref Q can be obtained. Our Pauli decomposition at 
threshold, eq. [T7] reveals the origin of some particularly interesting results. In particular, for later reference, we 
draw attention to the NR reduction of the amplitude of the topology with OgE within the ip (figs. 1(a) and 1(b), 
equation I13p. which involves the momentum q2 flowing through the photon-cc vertex and along the virtual fermion 
line connecting photon to gluon. The photo-produced cc pair is then clearly not in a simple ^Si state. Physically, 
contributions other than '^Si would vanish if the amplitude were proportional only to wavefunction at the origin; it is 
the spatial propagation away from the 'origin', associated with the propagator, that enables the non-zero amplitudes 
associated with these other configurations at threshold. 

First consider the factorizing contribution in eq. I17[ namely e^.ae^.Sj. The relative size of amplitudes arising 
from this term alone after combining with suitable Clebsch-Gordan coefficients in order to give helicity amplitudes 
for various x.i production in association with ip are as follows. 



V'(-)X2(++) = 1; V(0)X2(+) = l/\/2;^(+)x2(0) = 1/V6 (21) 
V'(0)xi(+) - l/V2;^(-t-)xi(0) = I/V2 (22) 
^(+)Xo(0) = 1/V3 (23) 

Thus for the unphysical case 5 = —1/2 where this term alone is present, one finds a result consistent with factorization 
as expected: 
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(^(to) : a(Vxi) : ^(V-Xo) = 5:3:1 (24) 

In the particular limit of threshold and 6 = 0, only the "non-factorizing" contribution £2-^1, = ' survives. 
This limit can be realised if the produced quarks have masses — > 0, as in the production of charmed mesons, D* , Dj. 

The expectation values of this term, after combining with the Clebsch-Gordan coefficients for Lx S J appropriate 
for XJ give, in the same normalisation as above 

^(+)x./(0) = [(J0|11; 1 - 1) + (J0|1 - 1; 11) - (J0|10; 10)] 

which vanishes for J = 1, 2 and is non-zero for J = where the constructive interference of the three terms gives the 
large value of a/3. 

For the case of e'^e~' i^ + X-i the amplitude is given by eq.[T9l The non-factorizing £2-££ = ^x'~^x term dominates 
by a factor of 4 relative to £i.£2 = • S^. While the factorizing term alone gave amplitude l/-\/3, the combination 
of the two terms now gives "^(1 ~ 12). Thus in leading order of NRQCD at threshold, we find the following relative 
rates: 

<e+e' ^ ^[X2 ■■ Xi;Xo]) = 5:3: 121. (25) 

This agrees with the result in ref 0; (to extract the OgE contribution, set F = in the amplitudes eqs. A3 in 
appendix A of ref Q)- 

The same combination of Clebsch-Gordan coefficients arises in the longitudinal amplitude for e^e~ —> ip{0)xjiO). 
Thus the origin of the large scalar amplitude at threshold is due to a S,,- • L^^, transition operator between the initial 
state and the L = 1 xj state which vanishes for all but J = 0. 

Note that the "0X1 ^-nd V'X2 amplitudes at threshold factorize. If the '0X2 amplitudes can be isolated, there is 
an interesting test concerning the non-zero V'(~)X2(++)- In ref [3| (table III and eqs. 36 — 41), it was noted for a 
^Si initial state that factorization implies that the 0'(— )X2(++) amplitude vanishes in S-wave, and that the D-wave 
amplitudes destructively cancel. In OgE however, there is more than simply "^Si in the initial state with the result 
that the S-wave is non-zero and as noted above we find 

a[0(-)X2(++)] : a[0(O)x2(+)] : a[V^(+)x2(0)] = 1 : 1/V2 : 1/V6 (26) 

in accord with S-wave dominance. 

Thus the 0'(— )X2(H — 1^) amplitude, far from vanishing, is predicted to be the dominant helicity state for ipX2- This 
provides an interesting test for the presence of non '^Si contributions in e+e^. 



C. One gluon exchange amplitude for arbitrary momentum 

Eq. [TH]is the amplitude for arbitrary r. From these amplitudes ref [t^ obtains rates for ip{Xi) + Xj{^2), where Ai 2 
are the helicities of the charmonium states, as a function of = IGni^/s. At the 10.6 GeV cm. energy of the data [g], 
where = 16m^/s = 0.28, ref f?\ finds for the OgE contribution to the cross sections 

a(0X2 : ^Xi : 0Xo) ^ 3 : 2 : 12. (27) 

The L and S dependence factorizes in three of the terms; the only term coupling L and S is the one already 
identified. Compared to the results at threshold, the relative sizes of (j{'ipX2) ■ o'('0Xi) have not changed much but 
there is an order of magnitude relative reduction of critpXc) at the higher energy. Nonetheless, it is the non-factorizing 
term that continues to dominate. Relatively large scalar meson production is predicted as a robust phenomenon at 
all energies. 



D. Charm pair production in NRQCD 

The above analysis can be applied to charm pair production, e+e^ ^ DjD^*^ or DjD^*\ Near thresh- 
old, factorization predicts 0| that the ratios of cross-sections from initial ^Si,^Di or vector hybrid to the states 
D*D2 : D*Di{^P) : D* Dq are as shown in table U 
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D*D2 


D*DiCP) 


D*Do 


'Si 





4 


3 




3 


5 





hybrid 


6 


3 


4 



TABLE L Relative sizes of transversely polarized decay amplitudes from OZI in charm pair production, 5^0. 

As in the cc +cc case we see here too that OZI factorization requires that the scalar production cannot be larger 
than both the axial and tensor channels. 

The interesting feature for charm pair production is that the OZI, or "flux-tube breaking" , mechanism is dynamically 
allowed, but that OgE can also be anticipated to be present. The relative sizes of these OgE contributions to the 
various charmed meson channels are calculated analogously to before, except that now we have S —f 0. Within the 
approximation that the large photon produces the cc, and the gluon then produces qq, the matrix element for 
D*Dj is given by eq. 1201 In practice there will also be some contribution where 7* — > qq followed by g — > cc. This is 
generally expected to be small because of unfavorable energetics. 

It is interesting to note that the £*.a,£j^.£2 structure of eq.l^Olis the same structure that, in eq.fTHl was responsible 
for the large 4'Xo amplitude in double charmonium production. Within the above approximation for D* Dj production, 
we see that the D*Dq channel is the only one driven by OgE, a maximal violation of factorization. The DDi{^Pi) 
channel is also non-zero. 

The double charmonium production is in practice the "worst case" for the scalar dominance. Yet it is clearly visible 
in the data, and our analysis at threshold shows why. For charmed meson production, where g qq with light 
flavours and S 0, the effect will necessarily be bigger, effectively infinite for reasonable parametrisation of mq/rric. 

The effective absence of OgE contributions for e+e^ — > D*Di and D*D2, which becomes exact in the limit S ^ 0, 
implies that the factorization selection rules should be particularly robust here. Thus we predict that at threshold, 
the helicity amplitude Z)*(±1)£'2(t2) = 0, in contrast to the case of '0(±1)x2(t2) where it is the largest of all the 
■0X2 amplitudes. 

The OgE selection of D*Dq and DDi in the threshold region may play some role in generating the enhancement 
seen as y(4260) (where charmed mesons in a relative S wave might rescatter to form ^pTnr). We shall return to 
the 7* qq process in more detail elsewhere. 



IV. CONCLUSIONS 

The intrigue of e+e" ^ tjj + ^ can hardly be overstated. There is a clear spectrum of C = H- charmonium states, 
whose pattern is not yet explained and in addition an enigmatic structure around 3940 MeV. As this is in the mass 
region where C = + charmonium hybrids could lurk [l^ , any theoretical modelling of this requires first understanding 
the population of the other bumps. The XJ production in particular needs to be understood; it seems empirically 
dominated by xoi a-t least when a mass fit is made to the cross section; no analysis has been made. We have drawn 
attention to the fact that the dominance of xo contrasts radically with normal OZI expectation where xo would be 
relatively small. We demonstrated that OgE gives large scalar, confirming Braaten and Lee 0], and have exposed the 
origin and significance of this result which had not previously been recognised as a sharp test of dynamics. 

Thus if dominance of "0X0 is confirmed over a range of q^ away from threshold, this would support OgE as the 
dominant decay mechanism. Conversely, if data near threshold confirm a[y(— )T(-|--|-)] —^ 0, this would signal factor- 
ization being dominant. In any event, we anticipate that the relative populations and helicity structures of 4'XJ will 
vary with q^. We recommend that this be investigated in e+e^ annihilation at super-B factories by means of ISR to 
access a range of energies. In particular experiment should attempt to measure the spin dependence of e~'"e~ ipX2 
as a function of q^ and compare with the analogous amplitudes in e+e^ 1^/2 or e^e~ — > p/2. 

When applied to charm pairs we find violation of OZI rules here too: OgE selects the D*Dq and DDi channels 
near threshold. The charm pair arena is interesting as it potentially enables us to test the relative role of OgE versus 
OZI dynamics. It also shows a possible mechanism for generating enhancements in these S'-wave production channels 
which, by constituent rearrangement, may produce these structures, such as seen in e+e" ipiriT at 4260 MeV [12|]. 

Our approach has demonstrated the breakdown of factorization in OgE. Ackleh, et. al. implicitly noted this in 
ref [ll| at least for on-shell constituents with a definite "^^^^Lj initial state. Our work goes beyond this by exposing 
the important role propagator effects can play; in particular we have found this to be essential in matching to the 
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NRQCD work of ref 0, at least at threshold. 

As a result, we have exposed the origin of the large amplitude for tpxa relative to V'Xi.2 which manifestly violates 
factorization and for which there are preliminary hints in data. We urge experiment to use spin analysis to confront 
this and other tests in both charm pair and light hadron production in order to sharpen understanding of the 
relative importance of OZI and non-factorizing dynamics. 
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